Draft version February 5, 2008 

Preprint typeset using I£T^X style cmulateapj v. 10/09/06 



DETONATING FAILED DEFLAGRATION MODEL OF THERMONUCLEAR SUPERNOVAE II. COMPARISON 

TO OBSERVATIONS 

Daniel Kasen 1,2 , Tomasz Plewa 3,4 

Draft version February 5, 2008 



o 
o 

(N 

o 

Q 



> 

00 
On 

(N 

o 



Of 
I 1 

o 



X 



ABSTRACT 

We develop and demonstrate the methodology of testing multi-dimensional supernova models 
against observations by studying the properties of one example of the detonation from failed deflagra- 
tion (DFD) explosion model of thermonuclear supernovae. Using time-dependent multi-dimensional 
radiative transfer calculations, we generate the synthetic broadband optical light curves, near-infrared 
light curves, color evolution curves, full spectral time-series, and spectropolarization of the model, 
as seen from various viewing angles. All model observables are critically evaluated against examples 
of well-observed, standard Type la supernovae (SNe la). We explore the consequences of the intrin- 
sic model asphericity by studying the dependence of the model emission on viewing angle, and by 
quantifying the resulting dispersion in (and internal correlations between) various model observables. 
These statistical properties of the model are also evaluated against those of the available observational 
sample of SNe la. On the whole, the DFD model shows good agreement with a broad range of SN la 
observations. Certain deficiencies arc also apparent, and point to further developments within the 
basic theoretical framework. We also identify several intriguing orientation effects in the model which 
suggest ways in which the asphericity of SNe la may contribute to their photometric and spectroscopic 
diversity and, conversely, how the relative homogeneity of SNe la constrains the degree of asymmetry 
allowable in the models. The comprehensive methodology adopted in this work proves an essential 
component of developing and validating theoretical supernova models, and helps motivate and clearly 
define future directions in both the modeling and the observation of SNe la. 

Subject headings: hydrodynamics — radiative transfer — spectropolarimetry — supernovae:general 



1. INTRODUCTION 

The critical assessment of theories of thermonuclear, or 
Type la, supernova (SN) explosions must ultimately in- 
volve not only the modeling of the complex 3-dimensional 
physics characterizing the explosion process itself, but 
also a comprehensive treatment of the multi-dimensional 
radiative transfer physics determining the emission from 
the ejected stellar material over the months and years fol- 
lowing explosion, thereby allowing the model to be tested 
by direct comparison to a broad range of supernova ob- 
servations, such as light curves, spectra and spectropo- 
larization, the aim being to reproduce the essential fea- 
tures of the observational sample, within the framework 
of a certain theoretical paradigm, and with the minimum 
number of free parameters and assumptions. 

Most theories of Type la supernova (SN la) explo- 
sions begin with the assumption that these events rep- 
resent the thermonuclear disruption of carbon-oxygen 
white dwarfs, typically near the Chandrasekhar limit. 
Despite several uncertainties concerning the structure of 
the progenitor star and the precise conditions leading 
to its ignition, complex multi-physics codes have been 
used to simulate the hydrodynamics and turbulent nu- 
clear combustion leading to the unbinding of the white 
dwarf. Nuclear burning typically lasts for up to 5 sec af- 
ter ignition; the stellar material ejected in the explosion 
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reaches the free-expansion phase ~ 1 minute later. At 
this time, the model output (i.e., the density, velocity, 
and compositional structure of the SN ejecta) may be 
fed into detailed radiative transfer codes for calculation 
of the synthetic model spectra, light curves and polariza- 
tion. 

The diagnostic value of such transfer calculations in 
assessing the model predictions is impressive. Light 
curve observations constrain the explosion energy, the 
total ejected mass, and the amount and distribution of 
the synthesized radioactive products (in particular 56 Ni) 
which power the SN light curve. Spectroscopic obser- 
vations constrain the ejecta velocity scale, thermal state 
and chemical stratification. Spectropolarization observa- 
tions constrain asymmetries and inhomogeneities in the 
ejecta density and compositional structure. Given the 
extreme complexity of the underlying phenomenon, no 
single measure can be used to establish the viability of a 
SN explosion model, rather one must consider this broad 
set of model observables as a whole. 

Inevitably, the comparison of model and observation 
involves a number of important uncertainties, including 
the physical approximations and geometrical simplifica- 
tions adopted in the theoretical calculations, as well as 
errors in the observations of SNe, their distances and the 
degree of dust extinction. Nevertheless, radiative trans- 
fer studies of 1-dimensional (ID) SN la models have al- 
ready provided importa nt insights into the nature and 
theory of these objects. iHoflich fc Khokhlovj (|1996t ). for 
example, synthesized light curves for a wide range of ID 
theoretical models, including deflagrations, pulsed deto- 
nations and delayed-detonations, and compared the re- 
sults directly to observations. The spectroscopic prop- 
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er ties of the paramete rized ID deflagration model w7 
of iNomoto et al.l (119841) have been frequently evaluated 
against observations (Branch et alJll98aiHarkness)[il)M 
Jefferv et a"llll992t iNugent et al.lll997t iLentz et alJboOl : 
Baron et al.l [2006). In addition, the synthetic light 
curves and spectra of ID delayed-detonation models have 
been studied i n great detai l, both in the optical and 
near infrared (iHoflicbl [19951 : IHoflich et all [19951 l200l 
iWheeler et al.lll998l ). All of these studies, by specifying 
the particular successes and failures of the models, have 
helped define the basic features required of any viable 
theory of SNe la. To date, however, there does not exist 
a self-consistent, fully unparameterized SN la model that 
satisfies all the observable properties of the class. 

With the advance of computing power, more realis- 
tic explosion models of SNe la have shown that multi- 
dimensional effects like ly play an essential role in the 
explo sion process (e.g., (Khokhlovl Il994t iReinecke et al.l 
2002). Corresponding advances in the observations, 
in particular polarization measurements, have given di- 
rect evidence that the ejecta of SNe la are not entirely 
spher i cally symmetric ([Howell et al.l 120011 IWang et al.l 
120031 iLeonardet all 120051 : IChornock et al.ll2006h . The 
validation of SN models faces new and interesting 
challenges once multi-dimensional scenarios are consid- 
ered. Although one can always study the gross, angle- 
averaged prope rties of aspherical mod els using ID trans- 
fer codes (e.g.. iBlinnikov et al. 2006), ultimately multi- 
dimensional transfer calculations are needed to retain the 
full geometry of the problem. In that way we can study 
how the model properties depend on viewing angle, and 
quantify the resulting dispersion in (and internal cor- 
relations between) the various model observables. Such 
statistical properties of the models may then also be eval- 
uated against those of the available observational sample 
of SNe la. 

In this paper, we develop and demonstrate such 
a methodology for testing multi-dimensional super- 
nova models by studying the observable properties of 
one example of the (axisymmetric) deto nating-failed- 
deflagration (DFD) model of SNe la (jPlewal I2006L 
hereafter Paper I). Using the multi-dimensional time- 
dependent radiative transfer code SED0NA, we calculate 
the full time evolution (day 2 to day 80 after explosion) 
of the emergent model spectrum as seen from various 
viewing angles. Broadband UBVRIJHK light curves 
are then generated by convolving the spectral series with 
the appropriate filter transmission profiles. The intrinsic 
asymmetry of the model is studied by calculating syn- 
thetic spectropolarimetry and by examining the depen- 
dence of all model observables on orientation. The high 
spectral resolution of the calculations allows us to study 
individual absorption features in the synthetic spectra, 
and in particular the level of line polarization and the 
temporal evolution of line Doppler shifts. All model 
properties are critically contrasted with representative 
examples of normal SNe la. 

Our approach constitutes a comprehensive, end-to-end 
simulation of the supernova evolution - from the first 
moments of ignition, through the weeks and months fol- 
lowing the explosion - without resorting to any ad-hoc 
manipulation of model parameters. This methodology 
provides the most thorough means of evaluating the suc- 
cess and failures of the particular explosion paradigm, 



and offers important guidance for further development 
of the theory. 

2. DFD Y12 EXPLOSION MODEL 

The detonating failed deflagration (DFD) scenario 
of thermonuclear supernovae considers an off-center, 
mild ignition process in a degenerate Chandrasekhar 
mass carbon-oxygen white dwarf. In Paper I, we de- 
scribed detailed calculations of the explosion process as 
perform ed with the multi-p hysics hydrodynamics code 
FLASH ([Frvxell et al.ll200o( ). These calculations were all 
performed in two dimensions assuming axial symmetry. 
From the set of the models presented in that work, we 
have selected one example (Y12) for detailed analysis. 
Here we briefly review the salient features of the Y12 
model from the point of view of the calculation of model 
observables. 

In the Y12 model, the white dwarf is initially ignited 
within a small spherical region on the symmetry axis, 
50 km in size and offset 12.5 km from the center. The re- 
sulting deflagration takes the form of a buoyantly rising, 
high temperature bubble filled with freshly synthesized 
nucleosynthetic products. This early phase is illustrated 
in Figure [lji, in which the bubble surface (identified with 
the edge of the deflagration front) is marked by the dark 
contour line. 

As the rising bubble approaches the white dwarf sur- 
face, the outer layers of the star are deformed and lat- 
erally accelerated, leading to the production of a circu- 
lar surface wave. By t=1.75 s (Figure [Tp) that wave 
is seen to already have passed the stellar equator. By 
t=3.575 s (Figure Ht) the wave reconverges at the pole 
and begins radially compressing unprocessed stellar ma- 
terial. A jet-like flow develops in this confined region, 
radially expelling material from the system while at the 
same time forming a shock wave moving into the bulk of 
the star. The shock-wave is marked by the light blue line 
near (r, z) ks (0, —2.5 x 10 8 ) cm in Figure Ht- 

A short time later, the temperature in that shock- 
wave-dominated region exceeds the ignition tempera- 
ture for the carbon/oxygen mixture, triggering a deto- 
nation through the shock-to-detonation mechanism (Pa- 
per I ). Unlike in previous delayed-detonat ion models 
(e.g., IHoflich et alj|1995t iGamezo et~aH l2005). this tran- 
sition to detonation is not inserted by hand but arises 
naturally in the course of the model evolution. Over 
the next 0.3 s, the detonation wave will progressively 
consume the bulk of the star (Figures [TJi and (Uj) , and 
reprocess the material ejected in the initial deflagration 
phase. Burning ceases a few tenths of a second later. 
By t — 5 s (Figure QJ) , the supernova Shockwave has al- 
ready completely swept through the portion of the star 
pre-expanded by the deflagrating bubble [z > 0), while 
it still remains trapped in the denser portion of the star 
(z<0). 

We followed the further evolution of the Y12 model up 
to t = 100 s after ignition. At this time, the ratio of in- 
ternal to kinetic energy is found to be less than 1% at all 
points in the ejecta, indicating that the free expansion 
phase has been reached. We verified that the model ve- 
locity structure is homologous, that is, the local velocity 
is linearly related to the radius v(r,t) — r(t)/to where 
to = 100 s. Thus spatial coordinates can be transformed 
to velocity space using v(r,t) = 10~ 7 f(t) km s _1 . The 
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Fig. 1. — Y12 model of a thermonuclear supernova explosion. Shown is the temporal evolution of density in log scale at times 
t = 1.25,1.75,3.575,3.75,4.0,5.0 s from panel (a) to (f) respectively (note the change of spatial scale between panel (f) and the rest of 
the sequence). The position of the deflagrating material is marked with black contours, while Shockwaves are marked in light blue. The 
transition from deflagration to detonation is triggered by the Shockwave visible in panel (c) at (r, z) = (0, —2.5 X 10 8 ). For more details, 
see the text. 



total kinetic energy of the explosion is 1.36 x 10 51 ergs 
and the total mass of 56 Ni produced is 0.93 M@. 

The density structure of the model at this final time 
(Figure \2§ is nearly spherical, but with a slight "egg- 
shape" , being more compact on one side. As seen in 
Figure [3l the characteristic ejection velocities are higher 
on the side from which the bubble emerged (the ignition 
side) and lower on the side from where the detonation 
took place (the detonation side). This is a result of the 
early expansion that occurred preferentially on the igni- 
tion side due to the passage of the deflagration bubble. 

The compositional structure of the model, also illus- 
trated in Figure El is well stratified on the large scale, 
and nearly spherical in shape. Interestingly, the entire 
compositional structure is offset by ~ 3000 km s _1 rel- 
ative to the center of expansion. The innermost region 
of ejecta is rich in iron-group elements (specifically 56 Ni) 
and extends to 13,000 km s _1 on the ignition side and 
7000 km s _1 on the detonation side. Surrounding this is 
a shell rich in intermediate mass elements (IMEs), span- 



ning the velocity range 7000 - 10, 000 km s~ 1 on the det- 
onation side and 13, 000 — 18, 000 km s _1 on the ignition 
side. Both 28 Si and 40 Ca are present in this IME shell, as 
marked respectively by the black and blue contour lines 
in Figure [2l Nearly the entire star is consumed by the 
explosion, with very little carbon/oxygen remaining. 

We note that the detailed nucleosynthesis calculated 
in the Y12 model is only approximate (see Paper I). 
During the deflagration phase, nucleosynthesis was cal- 
culated using an approximate burner which tracked only 
the few most representative isotopes. The detonation 
phase employed a more extensive (but still approximate) 
13-element network. Post-processing of the model using 
a more complete nuclear network is therefore necessary to 
more accurately specify the detailed model abundances. 
In particular, some fraction of the 56 Ni produced in the 
detonation phase may in fact represent stable iron group 
elements, while some fraction of the silicon group pro- 
duced in the deflagration phase likely represents a distri- 
bution of various IMEs. 

The inner layers of ejecta (v < 15,000 km s _1 ) in 
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Fig. 2. — Final density structure of the Y12 DFD model of a 
thermonuclear supernova explosion in the homologous phase (t = 
100 s). The density is shown in log scale. The spatial coordinates 
can be converted to velocity space using v(r, t) = 10~ 7 r(t) km s — 1 . 
The distributions of 28 Si and 40 Ca are shown in red and black 
contours, respectively. The contour levels correspond to Y( 28 Si) = 
0.25 and Y( 40 Ca) = 0.01. The position of the inner IME-rich 
shell can roughly be identified with the distribution of 40 Ca. The 
material in the outer silicon-rich region was produced during the 
deflagration phase. The lack of calcium in the outer shell is the 
result of the simplified nuclear burning method used to calculate 
the deflagration. Note the relative smoothness of the inner shell 
compared to the inhomogeneity of the outer shell. 

the Y12 model are smooth and continuous, typical of 
the material processed solely by the detonation. This 
can be contrasted with the outermost layers of ejecta, 
which were perturbed initially by the surface wave cre- 
ated by the deflagration bubble. That turbulent-like flow 
involved mixing of the bubble material and unburned 
stellar matter, creating density and compositional inho- 
mogeneities with a typical size of about ten degrees. 

A characteristic feature of the model is the presence 
of two distinct shells rich in IMEs. The inner IME shell 
represents material synthesized by the detonation wave. 
The outer, high-velocity shell (v « 25000-35000 km s _1 ) 
is composed chiefly of the material created and expelled 
during the deflagration phase. This high-velocity shell 
resembles the pan cake structure describe d in the param- 
eterized models of iKasen fc Plewal (|2005f ). except that in 
the Y12 model the material has a larger covering factor 
and a more irregular structure. Note that only silicon, 
and not calcium, is present in the high-velocity shell, an 
artifact of the approximate burner used during the defla- 
gration phase. One expects that, had a more complete 
nuclear network been used, a significant amount of cal- 
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Fig. 3. — Density structure (solid lines, in log scale) of the Y12 
DFD model in the homologous phase (t = 100 s) shown versus 
velocity coordinates. The dashed lines show the best fit exponential 
law p(v) = po exp(— v/v e ). Left: Density along the z-axis (the 
symmetry axis) at r = 0. The ejecta is more extended for z > 
(the ignition side) and characterized by an exponential velocity 
scale Ve = 3290 km s — 1 . The ejecta is more compact for z < 
(the detonation side) with v e = 1900 km s . The peak density 
is also slightly offset from the expansion center. Right Density 
along the r-axis at z = 0. The density law is characterized by 
v e = 2575 km s^ 1 

cium would also be present. Rather than speculate as 
to that exact amount of calcium produced, we leave the 
compositi onal structure as is, re ferring the reader to the 
studies of IKasen fe Plewal (|2005l ). 

Before considering the observable properties of the Y12 
model, we mention a few limitations of the underly- 
ing explosion calculations. Although no tuning of the 
model parameters was required to obtain an explosion, 
the DFD mechanism itself is subject to uncertainties 
regarding the initial conditions, geometrical simplifica- 
tions, and various numerical limitations. In Paper I, 
for example, a successful shock to detonation transition 
(SDT) was not achieved for all ignition configurations 
considered. This may indicate a lack of robustness of 
the underlying explosion mechanism, or may simply be 
a consequence of the relatively low numerical resolution 
of the calculations. A low resolution limits the size of 
the shocked regions, enhances numerical diffusion, and as 
a result may underestimate the actual post-shock tem- 
peratures and densities. These factors may artificially 
decrease the probability of SDT. On the other hand, nu- 
merical diffusion is a known source of mixing that may 
lead to artificial preheating and spurious ignition of fuel 



ig and sr. 

(|Frvxell, Miiller, fe Arnettj 11989). Although the transi- 
tions to detonation reported in Paper I are due to shock 
compression in nominally unburned material, the results 
are undoubtedly sensitive to possible contamination due 
to mixing, which may affect the likelihood of SDT. Even 
if a SDT is not achieved, it is still possible for a detona- 
tion to occur through the Zeldovich gradient mechanism 
(|Khokhlov. Oram fc Wheeler! fl997h once the deflagrat- 
ing material perturbs and mixes with the surface layers 
of the star. 

For completeness, we also summarize some more fun- 
damental simplifications of the DFD explosion model. 
In particular, we took the progenitor to be a cold, 
static, non-rotating, unmagnetized, chemically homoge- 
neous white dwarf. Although such a model is consid- 
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ered "standard" in most modern SN la explosion stud- 
ies, theoretical work strongly suggests that most of these 
assumptions are in fact incorrect. The cores of mas- 
sive white dwarfs are th ought to be both convective 
and c hemically stratified (jArnettl 119691 : iHoflich fc Stein! 
2002) . Energy release from carbon burning likely occurs 
even before ignition, which may contribute to the pre- 
expansion of the progenitor and modify the energetics of 
the subsequent deflagration. Rotation of the white-dwarf 
can affect the structure of the convective core and may 
induce a large scale anisotropy of the explosion. More- 
over, rotational shear may to some extent change the 
evolution of the deflagration front and further modify 
the effective flame energetics, which may influence the 
white dwarf pre-expansion. All of these elements neces- 
sarily contribute to the previously mentioned numerical 
uncertainties of the DFD explosion model, and deserve 
careful consideration in the future. 

3. RADIATIVE TRANSFER CALCULATIONS 

The observable properties of the Y12 DFD model 
have been calculated using the multi-dime nsional time- 
depen dent radiative transfer code SEDDNA (|Kasen et al.l 
2006). Given a homologously expanding SN ejecta struc- 
ture, SEDONA calculates the full time series of emergent 
spectra at high wavelength resolution. Broadband light 
curves are then constructed by convolving the spectrum 
at each time with the appropriate filter transmission pro- 
file. SEDDNA includes a detailed treatment of gamma-ray 
transfer to determine the instantaneous energy deposi- 
tion rate from radioactive 56 Ni and 56 Co decay. Radia- 
tive heating and cooling rates are evaluated from Monte 
Carlo estimators, and the temperature structure of the 
ejecta determined by iterating the model to thermal equi- 
librium. 

Several significant approximations are made in SEDDNA, 
notably the assumption of local thermodynamic equilib- 
rium (LTE) in computing the level populations. In ad- 
dition, bound-bound line transitions are treated using 
the expansion opacity formalism (implying the Sobolev 
approximation) and an approximate two-level atom ap- 
proach to wavelength redistribution. Special care was 
taken for the calcium lines, which were assumed to have 
pu re scattering s ource functions for the reasons discussed 
in iKasenl (|2006l ). Note that SEDONA allows for a direct 
Monte Carlo treatment of line fluorescence, but due to 
computational c onstraints this func tionality is not ex- 
ploited here. See iKasen et al.l |2006) for a detailed code 
description and verification, as well as tests of the expan- 
sion opacity and two-level atom approximations. 

In the atmospheres of SNe la, the microscopic condi- 
tions assuring the establishment of LTE ionization and 
excitation are in fact not met. Nevertheless, a number of 
previous theoretical studies confirm the adequacy of LTE 
models in reproducing the qualitativ e spectral and pho- 
tometric properties of SNe la (e.g., Hoflich et all 119951 ; 
iBaron et all 119961 : iPinto k. Eastman! |2001| ). The LTE 
simplification can be expected to result in quantitative 
errors in the model light curves, especially at later times 
(~ 60 days after explosion) when non-LTE ionization ef- 
fects become increasingly significant. In addition, the 
model predictions in this paper are subject to the fur- 
ther inadequacies of the two-l evel atom framework and 
the available atomic line data (|Kasen et alj|2006t [Kasenl 



2006). We hope to reduce these uncertainties in future 
work by relaxing some of the simplifying assumptions 
underlying the radiative transfer calculations. 

The time-dependent light curve calculations in this pa- 
per were computed using the following numerical resolu- 
tion: Spatial: A 100 x 50 cylindrical grid extending to 
a maximum velocity of 30,000 km s -1 and with a cell 
length of 600 km s 1 ; Temporal: 100 time points span- 
ning the epochs day 2 to day 80 with logarithmic spacing 
Alogi = 0.175; Wavelength: Spanning the range 100- 
25000 A with constant resolution of 100 A. Extensive 
testing confirms the adequacy of this resolution for the 
problem at hand. A total of 10 9 photon packets were 
used in the Monte Carlo procedure. The emergent light 
curves and spectra were generated by collecting escaping 
packets into 50 separate viewing angle bins. 

The time-dependent SEDDNA light curve calculation 
provides the full spectral time series of the model, how- 
ever to further improve our synthetic spectrum results we 
performed further time-independent calculations at five 
select epochs using both higher resolution (200 x 100 spa- 
tial grid, 5 A wavelength resolution) and a large number 
(10 8 ) of packets. The emergent bolometric luminosity 
parameter for these static calculations was determined 
from the full light curve calculation. The static syn- 
thetic spectra achieve higher signal-to-noise, but other- 
wise show excellent agreement with the time-dependent 
spectra, except in small-scale features of size < 100 A. In 
addition, one static spectrum near maximum light (t exp = 
18 days) was run using 10 12 packets in order to study the 
< 1% polarization levels. 

4. COMPARISON TO OBSERVATIONS 

We present below the synthetic optical/near- infrared 
light curves, spectral time series, and spectropolariza- 
tion of the Y12 model, comparing and contrasting the 
results with observations of representative SNe la. The 
properties of observed SNe la vary somewhat from object 
to object. We make no attempt here to find an optimal 
correspondence with the Y12 model, rather we simply 
compare to the more well-known, well-studied objects. 

4.1. Light Curves and Color Evolution 

In Figure HI we show the synthetic UBVRIJHK light 
curves of the Y12 model as compared to SN 2001cl, 
a well observed and photometrically normal SNe la 
(|Krisciunas et al.l2003h . There is some uncertainty in the 
total dust extinction to SN 2001el, which is estimated be- 
tween A v — 0.35 — 0.77 depending on the method used. 
Here we adopt the values A v = 0.57 a nd R v — 2.88 
(jKrisciunas et al.l[2003HWang et al. 2003) and correct f or 
extinction using the dust laws of ICardelli et al.l (|1989f) . 

Given the intrinsic asymmetry of the Y12 model, we 
consider for now only the equatorial view (9 = 90°) 
of the ejecta, deferring discussion of the orientation ef- 
fects to From this angle, the model light curves 
reach a peak absolute magnitude of -19.57 in B, and 
-19.61 in V. In Figure |4] we have adopted a distance 
modulus to SN 2001el of 31.65 mag so as to align the 
observed peak magnit udes with those of the model. 
iKrisciunas et a l. (2003) deduce a significantly lower dis- 
tance of 31.29 ± 0.08 mag on the basis of the super- 
nova photometry. This suggests that the Y12 model is 



6 



Kasen & Plewa 



= -19 - 



CD 
T3 
3 



C 

o 



J3 
< 




■19 



-18 



cr -16 
o 



.o 

< 



-15 



-14 



■13 til 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



20 30 40 50 
Days Since Explosion 



70 




20 30 40 50 
Days Since Explosion 



70 



Fig. 4. — Synthetic UBVRIJHK light curves of the Y12 DFD model as seen from an equatorial view (9 = 90°, solid lines) compared 
to observations of SN 2001cl (filled circles, Krisciunas ct al. 2003). To align the models and observations, we adopt a distance modulus for 
SN 2001el of fi = 31.65, which is 0.35 mag greater than the best estimate. 



~ 0.35 mag too bright compared to SN 2001el - not sur- 
prising given the explosion produced over 0.9 Mq of 56 Ni. 
This discrepancy may be mitigated somewhat if detailed 
nucleosynthesis calculations reveal that some fraction of 
the 56 Ni in the Y12 model in fact corresponds to stable 
iron group elements. Note also that the DFD scenario 
encompasses a family of explosion models which produce 
greater or lesser amounts of 56 Ni, depending on the ini- 
tial ignition conditions (Paper I). 

The shape of the model optical light curves match qual- 
itatively those of SN 2001el, especially in the B-, V- and 
i?-bands. The rise time to B-band maximum is 18.8 days 
for the model, similar to that observed amo ng typical 
SNe la (|Riess et al.lll999aHConlev et al.ll2006t) . The sub- 
sequent decline from peak is also quite consistent with 
the observations of SN 200 lei and other typical SNe la, 
with a drop in S-band magnitude 15 days after peak 
of Ami5(i?) = 1.19 mag. However, given that brighter 
SNe la generally display broader light curves, this de- 
cline rate may be considered too fast in light of the high 
peak magnitude of the model. As we discuss in Sj5[ the 
decline rates of the model light curves in fact depend 



rather sensitively on the orientation. 

A distinct secondary maximum in the model I- and 
./-band light curves occurs about three weeks after B- 
maximum. A much weaker secondary maximum is also 
present in the i?-band. The same double-peaked mor- 
phology is seen in the SN 2001el observations, although 
the size and timing of the secondary maximum differ 
somewhat from the model. The properties of the sec- 
ondary maximum in fact vary substant ially from su- 
pernova to supernova (|Nobili et al.ll2005f ), and also de - 
pend sensitively on the model parameters (Kasen 2006). 
Brighter SNe la (i.e., those with larger masses of 56 Ni) 
typically display more prominent secondary maxima. 
The overly bright /-band secondary maximum of the Y12 
synthetic light curve indicates, in part, its overly large 
mass of 56 Ni compared to SN 2001el. 

The model H- and if -band light curves show a much 
poorer correspondence with the observations, with the 
secondary maximum being m uch too brigh t relative to 
the first one. As discussed bv iKasenl (|2006h . this sort of 
discrepancy is attributable, in large part, to the inade- 
quacy of the atomic line data used in the transfer calcula- 
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Fig. 5. — Temporal evolution of the colors of the Y12 DFD model 
as seen from an equatorial view (6 = 90 °, solid lines) compared to 
observations of SN 2001el (open circles, Krisciunas ct al. 2003). 

tions. The Y12 model, however, shows a relatively poorer 
match to SN 2001el wh en compared to other models con- 
sidered in lKasenl (|2006l ). likely due again to the relatively 
large 56 Ni mass. 

The optical color curves of the model are compared to 
those of SN 2001cl in Figure [5l In general, the curves 
show qualitative agreement up to about 40 days after 
i?-band maximum. At later times, the model's color get 
progressively redder, while the observations show the op- 
posite trend. It is at these late epochs, when LTE pre- 
dicts neutrality in the ejecta, that non-thermal ionization 
by the produc ts of radioactive decay becomes significant 
(|Swartdll99ll) . Thus our transfer calculations should be 
considered less reliable at these phases. 

4.2. Spectral Evolution 

The spectra of SNe la consist of broad P Cygni fea- 
tures superimposed on pseudo-blackbody continuum. By 
identifying individual absorption features and measuring 
their Doppler-shifts, we probe the composition and ve- 
locity of the absorbing ejecta material. As the SN mate- 
rial expands and geometrically dilutes, the photosphere 
recedes, revealing emission from progressively deeper lay- 
ers. The spectral time-series thus provides a scan of the 
multi-layered ejecta compositional structure. 

We choose to compare the synthetic spectra of the 
Y12 DFD model to observations of SN 1994D, a well- 
sampled SN la with normal spect ral characteristics 
(jPatat et al.lll996HMeikle et al.lll996h . The dust extinc- 
tion of SN 1994D is believed to b e very small, and we 
make no corrections for reddening. iBranch et alj (|2005f ) 
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Fig. 6. — Synthetic spectra at five different epochs of the Y12 
DFD model as seen from an equatorial view (8 = 90°, black 
lines) compared to those of SN 1994D (red lines, Patat ct al. 1996; 
IMeikle ct al. 1996). Labels give the day relative to B-band maxi- 



provide a comprehensive line identification and spectro- 
scopic analysis of SN 1994D. In addition, SN 1994D has 
been used to validate sev eral other 1-dimen sional explo- 
sion models such as W7 (iLentz et~aL 1200 1| ) and various 
delay ed-detonation models ( Hoflichl 19951 : iBaron et ahl 
2006), allowing cross-comparison of the Y12 model with 
other theoretical calculations. 

In Figure [6l we show the synthetic spectra of the Y12 
model (as seen from an equatorial view) at several differ- 
ent epochs - viz., day —7, 0, +7, +15, and +25 days with 
respect to B-band maximum. In general, we find the 
synthetic spectra to be in impressively good agreement 
with those of SN 1994D. The quality of fit is compara- 
ble to (and often better than) that obtained with clas- 
sic 1-dimensional models, such as W7 and the delayed- 
detonation models. 

At maximum light (t cxp = 18.8 days), the model re- 
produces all of the major spectral features identified in 
the observations. These are primarily lines from IMEs, 
in particular three Si II absorption features (near 4000 A, 
5850 A, and 6150 A), the S II "W" -feature (near 5200 A), 
and the Ca II H&K and IR-triplet features (near 3600 A 
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Fig. 7. — Variation with viewing angle of the Y12 DFD model 
synthetic spectrum at maximum light (t oxp = 18 days). The angle 
8 = 0° corresponds to the view from the ignition side (z > in 
Figure [2}- 

and 8200 A, respectively). The only feature not well fit 
at maximum light is the double absorption near 4600 A, 
due to blends of Mg II and Si III. Interestingly, this fea- 
ture is peculiar in SN 1994D, and most other normal 
SNe la more closely resemble the model (e.g., SN 1981B, 
iBranch et aflllQSffl . 

The blueshifts of the absorption minima in the max- 
imum light spectrum indicate that the model velocities 
are in broad agreement with the observations. In de- 
tail, however, the model absorptions are systematically 
bluer than the observations, most conspicuously in the 
Si II 6150 and Ca II IR triplet lines. We note that the 
blueshift of the model absorption changes with orienta- 
tion (see a nd also is observed t o vary from supernova 
to supernova (|Benetti et alj I2005T ). The velocity differ- 
ences noted in Figure [5] are therefore within the observed 
scatter of normal SNe la. 

A week after maximum light (t exp — 25 days), most 
of the same IME features are still present in the ob- 
served spectra, while a prominent Na I feature emerges at 
5600 A, blending with the smaller of the Si II features. 
The approximate alpha network used in the explosion 
model did not include sodium, therefore the synthetic 
spectrum does not reproduce the sodium feature. This is- 
sue will be addressed in future models by including more 
detailed nucleosynthesis. Among the other features, we 
note that the sulfur W-feature, well fit by the model at 
maximum light, appears too weak at this epoch. 

At later times, the SN photosphere recedes into the 
metal core of cjccta, and the spectra progressively be- 
come dominated by lines from iron group elements. The 
model spectra at 14 and 21 days after B-maximum well 
account for the observed features, with the exception of 
the Na I line. 

5. INTRINSIC MODEL ASYMMETRY AND DIVERSITY 

A major advance offered by our multi-dimensional ra- 
diative transfer calculations is the ability to study the 
dependence of the model observables on the viewing an- 
gle. In the DFD model, the fundamental orientation ef- 
fects arise from the variation of the characteristic ejecta 
expansion velocities along different directions. As dis- 
cussed in $2] (see also Paper I) this distortion reflects the 



Fig. 8. — Temporal evolution of the Si II 6150 absorption min- 
imum velocity as seen from three different viewing angles. The 
angle 9 = 0° corresponds to the view from the ignition side (z > 
in Figure [SJ. while 6 = 90° corresponds to an equatorial view and 
8 = 180° corresponds to the view from the detonation side [z < 0). 
The variation in the model ve locity evolution can be compared the 
observed diversity in SNe la (Benctti et al. 2005, Figure 1). 

inherent asymmetry of the initial conditions and subse- 
quent deflagration. The global asymmetry of the par- 
tially burned white dwarf is then imprinted on the dy- 
namical and compositional structure of the ejecta by the 
off-center detonation. 

5.1. Photometric and Spectroscopic Diversity 

The asymmetry in the Y12 ejecta structure leads to 
two significant orientation effects in the model spec- 
trum, both clearly visible in the maximum-light spec- 
trum shown in Figure [7] 

First, the blueshift of most spectral absorption features 
depends on the viewing angle, becoming larger as one 
looks closer to the ignition side. This is due to the gen- 
erally higher expansion velocities of the IME shell on that 
side of the ejecta. Figure [8] quantifies this Doppler shift 
dependence by plotting the time evolution of the veloc- 
ity of the prominent Si II absorption near 6150 A, as 
seen from several viewing angles. In this figure we have 
measured the velocity directly from the Doppler shift of 
the Si II absorption minimum for each of the synthetic 
spectra in our daily spectral time-series. Such a method 
corresponds exactly to what is done in the observational 
studies. By contrast, most other time-dependent SN ra- 
diative transfer calculations determine only the velocity 
of the SN photosphere, which may or may not coincide 
with that of the Si II line. 

Near maximum light, the Si II velocity varies from 
9,000 to 13,500 kms -1 , or by nearly 40%, depending 
on the viewing angle. The velocity is largest and has the 
fastest rate of decrease with time when observed near the 
ignition side (9 — 0°). For viewing angles away from the 
ignition side, the velocity of the Si II feature is lower and 
almost constant with time. This va riety of velocity grad i- 
ents can be compared to Figure 1 of lBenetti et~aTT (2005). 
which demonstrates a strikingly similar diversity among 
observed SNe la. The Y12 model prediction is there- 
fore, at first glance, in reasonable agreement with the 
statistical behavior of the observat ional sample. On the 
other hand, the cluster analysis of lBenetti et alj (|2005f ) 
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suggests that observed SNe la can be divided into two 
distinct classes of "high" and "low" velocity gradient. In 
the Y12 model, by contrast, the orientation effect leads 
to a range of events intermediate the two classes. This 
continuous behavior ap pears, at first g lance , inconsistent 
with the dichotomy of iBenetti et al.l (|2005h . although a 
detailed statistical comparison would be needed to truly 
confirm the discrepancy. 

The second consequence of the model asymmetry on 
the maximum light spectrum is that the flux level at 
shorter wavelength shows a strong variation with view- 
ing angle. In the B-band (4000 A < A < 5000 A) this 
variation is only of order 10%, however it increases to 
a factor of ~ 4 in the [/-band (3200 A < A < 4000 A) 
and is even larger in the ultraviolet. The effect is a re- 
sult of the aspherical distribution of iron group elements 
in the ejecta. On the ignition side (9 = 0°) iron group 
species extend to a velocity v w 13, 000 km s~ , whereas 
on the detonation side (9 — 180°) these elements are re- 
stricted to velocities v < 7,000 km s _ . It is primarily 
the line blanketing from iron group lines that determines 
the flux level at shorter wavelengths. When viewed at 
9 = 180°, the relative lack of iron group elements at 
higher velocities leads to significant reduction in the de- 
gree of line blanketing, and hence a spectrum that is 
brighter at bluer wavelengths. 

At two weeks after maximum light, the same sort of 
orientation effects are visible in the synthetic spectra 
(Figure [9]). By this epoch, the layers of ejecta rich in 
iron-group elements have begun to recombine from dou- 
bly to singly ionized, and features from Fe II/Co II lines 
dominate the spectrum. Again, because of the aspheri- 
cal distribution of iron group elements, the flux level at 
shorter wavelengths shows a strong dependence on the 
viewing angle. However, at this time the asymmetry af- 
fects wavelengths even longer than at maximum light, 
as Fe II/Co II lines have become prominent in the op- 
tical part of the spectrum. The £>-band portion of the 
spectrum now shows flux variations by a factor of nearly 
two. This line-blanketing effect leads to a strong depen- 
dence of the model B — V color evolution on orientation, 
especially after maximum light, as seen in Figure fTTJl 

Clearly this orientation effect has important implica- 
tions for the B-band light curves of the model, and con- 
sequently the SN la width-luminosity relationship. The 
development of Fe II/Co II line blanketing largely de- 
termines the .B-band de cline rate Ami^{B) of SNe la 
(|Kasen fe Wooslevll2006[ ). As seen in Figure fTTl the as- 
pherical distribution of iron group elements in the Y12 
model thus leads to significant variations in Ami5(B) 
with viewing angle. The decline rate changes from 
Am 15 (B)= 0.96 mag at 9 = 0° to Ami 5 (B)= 1.36 mag 
at 9 = 180°. This is a good fraction of the total range of 
decline rates noted in the observational sample of SNe la. 

The peak .B-band magnitude of the Y12 model also 
varies with viewing angle, although to a lesser ex- 
tent, from M B = -19.43 mag at 9 = 0° to M B = 
— 19.56 mag from the opposite side. The correlation be- 
tween Anii5(B) and Mb thus has the same sign as the 
observed width-luminosity relation, in that broader SNe 
are generally brighter. However, as seen in Figure 1111 
(right panel), the correlati on is too weak as co mpared 
to the observed relation of iPhillips et al.l (|1999f) . Thus 
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Fig. 9. — Variation with viewing angle of the Y12 DFD model 
synthetic spectrum two weeks past maximum light (t escp = 32 days). 
The angle 9 = 0° corresponds to the view from the ignition side 
(z > in Figure [2} Note that the strong flux variation has shifted 
redward when compared to that at maximum light (Figure 0. 
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Fig. 10. — Variation with viewing angle of the Y12 DFD model 
synthetic B — V color with time. The angle = 0° corresponds to 
the view from the ignition side (z > in Figure [2}- 

the asymmetry of the Y12 model could at most be con- 
sidered as a potential source of dispersion in the width- 
luminosity relation. The model intrinsic dispersion is 
cr w 0.13 mag, roughly equa l to the level noted in obser- 
vations (jHamuy et al.l fl996) and leaving little room for 
a contribution from observational or other errors. 

As also seen in Figure [11] (right panel), the light curves 
of the Y12 model appear to be systematically brighter 
than the observations, given their range of Ami5(£?). 
The shade d band in the figure is the width-luminosity 
relation of IPhillips et alj ([1999) with an assumed cali- 
bration of Mb — —19.3 mag for Amis(B) = 1.1 mag 
and with a dispersion a = 0.15 mag. The model points 
lie roughly 0.25 mag above this relation. The signifi- 
cance of this discrepancy is debatable, as the absolute 
calibration of the empirical width-luminosity relation is 
subject to uncertainties in the distances to SNe la, while 
the model light curves may possess systematic errors of 
order 0.1-0.2 mag due to approximations made in the 
radiative transfer calculations. 

The variation in the £?-band magnitude with viewing 
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Fig. 11. — Left: Variation of the Y12 model B-band light curve with viewing angle. 8 = corresponds to the view from the ignition side 
(z > in Figure[2]l. Right: Relationship between the model peak B-band magnitude Mb and decline rate A mi^(B) as seen from d ifferent 
viewing angles. Color coding is the same as in the left panel. The shaded region is the observed relation of Phillips ct al. (1999) with a 
calibration Mg = —19.3 mag at Amis(B) = 1.1 mag and spread of a = 0.15 mag. The dotted lines show the same relationship shifted 
upward by 0.25 mag. 



angle becomes even more dramatic at later times, and 
exceeds 1.5 mag at 40 days after B- maximum. This is 
larger than th e ~ 1 mag dispersio n seen in the observa- 
tional sample (|Riess et al.lll999"H) . Note, however, that 
the increased uncertainties in the radiative transfer at 
later times (as noted in i j4.1[) may contribute to the ex- 
treme spread in model B-magnitudes at these epochs. 

5.2. Polarization 

Polarization measurements provide the most direct 
probe of the asphericity of the SN ejecta. Light is po- 
larized by electron scattering in the hot, ionized SN 
envelope. Detection of a finite net polarization indi- 
cates a preferred direction in the scaterring medium, 
and hence a departure from spherical symmetry. Spec- 
tropolarimetry is a challenging endeavor, as the light 
from SNe is dim and quickly fading, and the observed 
polarization levels are typically small, P < 1%. Nev- 
ertheless, intrinsic polarization has now been clearly 
detected in sev eral SNe la, both in the lines and 
the continuum dHowell et all l200il [Wang et ail 120031 : 



iLeonard et all l2005t IChornock et at] 120061) . The sam- 
ple size of objects with high signal-to-noise spectropo- 
larimetry is unfortunately still rather small, and the ob- 
servations are also subject to uncertainties in the degree 
of interstellar polarization caused by scattering by dust 
grains along the line of sight. 

The continuum polarization level (i.e., the degree of 
polarization at wavelengths where electron scattering is 



the dominate opacity) provides the most basic measure 
of global asymmetry in the bulk SN ejecta. For the Y12 
model, the continuum polarization is in fact very small 
during the epochs around and after maximum light. This 
reflects the fact that the density contours of the ejecta 
are not grossly aspherical. Moreover, the primary mode 
of the model's "egg-shaped" density structure is on a 
180° scale, whereas the continuum polarization measures 
differences oriented 90° apart. In particular, while the 
polarization of light scattered from the ignition side of 
the ejecta is enhanced compared to a spherical model, 
this is nearly counter-balanced by the relative decrease 
in the polarization of light scattered from the compact 
detonation side 180° away. The polarized flux from the 
polar regions of the Y12 model then nearly equals the 
orthogonally oriented polarized flux from the equatorial 
regions. We therefore find that at maximum light, the 
continuum polarization of the Y12 model is small, P < 
0.1%, from all viewing angles. This is smaller than the 
continuum polarization of P ~ 0.3% typica lly observed 
in SNe la (e.g., SN 2001el. lWang et al.ll200l . Given the 
signal-to-noise of most spectropolarimetry observations, 
a polarization of P < 0.1% is usually indistinguishable 
from zero. 

Although the continuum polarization of the Y12 model 
is small at maximum light, the polarization over lino fea- 
tures can be substantial. In Figure [T2l we show the syn- 
thetic model polarization spectrum at maximum light, as 
seen from several viewing angles. Although nearly 10 12 



Observational Properties of DFD 



11 



e = 10° 




4000 5000 6000 7000 8000 9000 

Wavelength (Angstroms) 



Fig. 12. — Synthetic polarization spectrum of the Y12 DFD 
model near maximum light (t CX p= 18 days) as seen from several 
viewing angles 9. Small scale fluctuations represent the level of 
Monte Carlo noise in the calculation. 

photon packets were used in this Monte Carlo calcula- 
tion, the synthetic polarization spectra still show random 
noise at a level ~ 0.1%. Nevertheless, the broad line 
features are easily discernible among the high-frequency 
fluctuations. The most prominent polarization features 
are due to the Si II absorption feature near 6150 A, the 
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Fig. 13. — Velocity (open circles, in 10 4 km s — x ) and absolute 
value of the polarization (closed circles, in percent) measured at 
the line minimum of the Si II 6150 feature of the Y12 model at 
maximum light. The size of the error bars (0.1%) roughly repre- 
sents the level of Monte Carlo noise in the polarization calculations. 
The abscissa is x = (1 — cos 9)/2 where 9 is the viewing angle. The 
probability of observing the model within the range x = (3:1,3:2) 
is then easily figured as p = X2 — xi . The polarization cancels by 
symmetry at 9 = 0° and 180° and so has been set to zero at these 
points. 

Ca II IR triplet near 8200 A, and the blend of Fe II lines 
near 4900 A. The polarization over the Si II absorption 
feature reaches nearly 2% from certain viewing angles. 

The line polarization of the Y12 model reflects the 
global asymmetry of the ejecta compositional structure, 
which results in an asymmetrical distribution of the line 
opacity. From certain viewing angles, the opacity from 
a strong line partially eclipses the underlying electron 
scattering photosphere, leading to an incomplete can- 
cellation of the polariza tion over the wave lengths of the 
line absorption feature (|Kasen et~all[2003h . In the Y12 
model, the most important asymmetry is that the dis- 
tribution of IME-rich material extends to a higher veloc- 
ity on the ignition side of the ejecta. In some respects, 
this asymmetr y resembles that of the "pancake" geom- 
etry studied in iKasen &: Plewal (|2005f h although in the 
present case the excess of IME-rich material abuts the 
primary IME shell, rather than being detached from it. 
The dependence of the line polari zation on viewin g angle 
thus resembles that explained in IKasen fc Plewal (2005) 
- in particular, the line polarization is negative for small 
viewing angles (9 < 45°) but flips sign for intermediate 
viewing angles (9 w 90°). 

Although we have no direct way to determine the view- 
ing angle itself from SN observations, what we can study 
is the correlations between different angle dependent ob- 
servables. In particular, it is useful to study the relation 
between the velocity of a line and its polarization, where 
both quantities are measured at the minimum of the flux 
absorption feature. In Figure[T31 we plot these values for 
the Si II 6150 feature of the Y12 model, as derived from 
the maximum light polarization spectrum. In the figure, 
we chose the abscissa to be x = (1 — cos#)/2, where 9 
is the viewing angle. The probability of observing the 
model within the range x = (x\, X2) is then simply given 
by p = x 2 - x\. 

As discussed already ( q5.1[) , for viewing angles towards 
the ignition side (x < 0.1 or 9 < 35°) the Si II velocity is 
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unusually high (v w 14,000 km s -1 ). From these angles 
the Si II polarization is also always large (P ~ 1 — 2%). 
A similar trend appears in the observations. Among the 
limited number of SNe la with spectropolarimetry, the 
Si II polarization is found to be largest in the two objects 
with the highest Si II absorption ve locities (SN 2002bf 
and SN 2004dt. ILeonard et al.ll200l . The Y12 model 
predicts such events should be observed roughly 10% of 
the time. 

From viewing angles away from the ignition side, the 
Si II polarization of the Y12 model may be either small 
or large. For 45° > 9 > 78° (x = 0.15 - 0.4), the 
model generally reproduces the low polarization lev- 
els and moderate line velocities charact eristic of typical 
events like SN 200 le i (|Wang et al.lf2003h and SN 2003du 
(|Leonard et al.ll2005[ ). However, the Y12 model also pre- 
dicts that some SNe la should exhibit high Si II polar- 
ization accompanied by low or moderate Si II line veloc- 
ity. In particular, for x = (0.45,0.70) and x = (0.9, 1.0) 
the polarization is P < 1.0% while the velocity is only 
v w 9,000 - 10,000 km s -1 . According to the model, 
such objects should be observed roughly 35% of the time, 
however no SN la with this characteristic has yet been 
found. 

On the whole, the Y12 model predicts that the Si II line 
polarization will be greater than 1% from almost half of 
all possible viewing angles, which is generally higher than 
the average level of line polarization noted in the current 
(albeit limited) observational sample of SNe la. This 
suggests that the degree of compositional asymmetry in 
the model ejecta structure may be too extreme. This 
same conclusion had already been indicated by the large 
dispersion in the B-band decline rates ( £)5.1|) . 

In addition to a global asymmetry, the IME layer of 
the Y12 model possesses numerous compositional inho- 
mogeneities ("clumps") of typical size ten degrees. These 
were the result of the original turbulent-like flow of the 
deflagration material across the white dwarf surface. The 
signatures of these clumps, however, are not readily dis- 
cernible in the polarization spectra of Figure 1121 In this 
2D calculation, the global asymmetry of the IME distri- 
bution dominates the polarization signal. In a fully 3D 
calculation, on the other hand, the IME clumps would 
lead to some deviation of the line polarization from a 
single axis of symmetry. This spectropolarimetric sig- 
nature has been used to infer c lumpiness in the IME 
layers of some observed SNe la ()Wang et al.ll2004l ). In 
this sense, the DFD scenario is notable in its ability to 
produce a clumpy outer layer of IME while retaining a 
rather smooth distribution in the inner ejecta. 

The Y12 model also possesses clumps of IME at very 
high velocity (v = 25, 000 - 35, 000 km s" 1 ) representing 
material produced and expelled during the deflagration 
phase. Assuming these clumps were to possess a moder- 
ate abundance of calcium, one would expect to observe 
high velocity, highly polarized Ca II IR triplet absorp- 
tion features in th e maximum- ligh t spect rum, similar to 
those discussed in lKasen fc Plewal (j2005) . However, be- 
cause the approximate nucleosynthesis employed during 
the deflagration phase of the explosion calculation did 
not include calcium, the high-velocity material does not 
in fact influence our current synthetic flux or polarization 
spectra. Given the low densities in these layers, the lines 



from all other elements besides calcium are too weak to 
be seen. 

More detailed and conclusive studies of the ejecta ge- 
ometry of SNe la will be possible once the database of 
spectropolarimetry observations is expanded, such that 
for each subclass of SNe la we have enough events to 
sample all viewing angles in a uniform way (e.g., always 
at exactly the same epoch). Figure [T3l is just one demon- 
stration of how we can use such a data set to powerfully 
constrain multi-dimensional models of SNe la. 

6. SUMMARY AND CONCLUSIONS 

We presented the model observable properties of one 
example of the DFD explosion scenario of thermonuclear 
supernovae. The explosion model involved the off-center 
ignition of Chandrasekhar-mass white dwarf which later 
underwent a transition from deflagration to detonation 
and produced an energetic explosion. The model was ex- 
panded well into the homologous phase, and then used 
as input to time-dependent multi-dimensional radiative 
transfer calculations. The transfer calculations supplied 
the emergent broadband light curves, color evolution, 
spectral time-series, and polarization of the model, all 
of which were critically evaluated against examples of 
well-observed, standard SNe la. 

The calculations presented here were intended pri- 
marily to develop and apply a methodology of validat- 
ing multi-dimensional SN models against observations. 
They represent the first extensive comparison of a self- 
consistent multi-dimensional SN la model against obser- 
vations, and comprise one of the most comprehensive 
evaluations of any SN model to date. On the whole, 
the properties of the Y12 DFD model appear fairly con- 
sistent with observations, although certain discrepancies 
are also apparent. We summarize the primary successes 
and failures of the model as follows: 

• The shape of the model broadband light curves 
resembles observations of the normal Type la 
SN 2001el. However, given the model's large 56 Ni 
mass (0.93 M Q ), the synthetic light curves are gen- 
erally brighter than the typical SNe la, by about 
0.35 mag. The model light curves reach a peak ab- 
solute magnitude of -19.57 in B, and -19.61 in V, 
as viewed from the equatorial direction (6 — 90°). 

• The i?-band decline rate of the model (Amis(i?) = 
1.19 mag, as viewed from the equatorial direction) 
is typical of an average SN la, but too fast when 
compared to observed SNe la of equivalently high 
peak brightness. 

• The R-, I-, and J-band light curves of the model 
exhibit a secondary maximum which resembles 
qualitatively that seen in observations. The timing 
and prominence of the model secondary maximum 
differs somewhat from that seen in SN 2001el, likely 
due to the model's relatively large 56 Ni mass. 

• The B — V and V — R color curves of the model 
(from 9 = 90°) fit those of SN 2001el to within 
~ 0.1 mag up to day 30 after maximum light. 
The U- and /-band light curves more strongly de- 
viate from the observations, leading to a slightly 
poorer, but still reasonable match between the 
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model U — B, B — I, V — I, and R — I color curves 
and observations. 

Inspection of the model spectra (from an equato- 
rial view) each week from -7 to +25 days relative 
_B-maximum shows good agreement with observa- 
tions of the spectroscopically normal SN 1994D. 
The quality of fit is comparable to that of stan- 
dard ID explosion models such as W7 and the ID 
delayed-detonations . 

The asymmetry of the model ejecta structure leads 
to intriguing orientation effects. The blueshifts of 
individual absorption features vary with viewing 
angle by more than 40% at maximum light. More- 
over, depending on the viewing angle, the evolution 
of the expansion velocities (as measured from the 
minimum of the Si II 6150 absorption) accommo- 
date the entire range of obse rved behaviors, from 
high to low velocity gradient (jBenetti et al.l l2005). 

Due to the asymmetrical velocity distribution of 
iron group elements, first the ultraviolet and later 
the blue part of the spectrum show substantial 
variation (by a factor of several) in flux absorp- 
tion along different viewing angles. This leads 
to a significant dependence of the B-band decline 
rate on orientation, with ATO15 varying between 
0.96 — 1.36 mag. The model shows a correlation 
of Ami5(B) with the -B-band peak magnitude, but 
the amplitude is twice weaker than the observed 
Phillips relation. Thus in the present model, the 
orientation effects can be considered as a possible 
source of dispersion in the underlying relation, at 
the level 0.12 mag. 

The intrinsic mild asymmetry of the model den- 
sity structure, leads to a continuum polarization 
at maximum light of P < 0.1%, somewhat lower 
than that typically observed. This may indicate 
that other effects not included in the present model 
(e.g., rotation) may be also required to reproduce 
the large scale ejecta asymmetry. 

The model shows line polarization up to a level P ~ 
2% in the Si II "6150" and Ca II IR-triplet features, 
due to the aspherical distribution of IMEs. The line 
polarization is large from inclinations where the 
line velocities are highest, resemb l ing th e observed 
behavior noted in iLeonard et alj (|2005f ). However 
significant line polarization (P ~ 1 — 2%) can also 
be seen from some viewing angles where the line 
velocities are relatively low, a characteristic that 
as of yet has no correspondence in observations. 

In the model, the material synthesized during the 
initial deflagration phase is transformed into a high 
velocity {v ~ 25,000 - 35,000 km s _1 ) IME-rich 
outer shell. This shell resembles the structure 
considered in Ka sen fc Plewal (|2005l ). but in the 
present case has a larger covering factor and a more 
irregular structure. Although not considered di- 
rectly in this paper, this material offers one pos- 
sibility to explain the detached, high-velocity ab- 
sorption featur es commonly observ ed in early time 
SN la spectra (|Mazzali et al.ll2005D . 



The comprehensive methodology applied in this paper 
helped to appreciate certain deficiencies in the theoretical 
calculations. The failure to reproduce the Na I features in 
the post-maximum spectra follows from the use of an ap- 
proximate alpha-network nucleosynthesis scheme lacking 
sodium. Similarly, the lack of calcium in the outer high- 
velocity IME shell is a direct consequence of the crude nu- 
clear burning scheme used during the deflagration phase. 
The wrong trends noted at later time (t eKp > 60 days) 
model colors likely indicate a breakdown of the assump- 
tions adopted in the radiative transfer calculations, in 
particular the neglect of non-thermal ionization effects. 
Meanwhile, the poor correspondence in the model H- 
and -ftf-band light curves with observations reflects, in 
large part, the inadequacy of the available atomic line 
list data. 

The current Y12 model is only one example of the fam- 
ily of DFD explosions (Paper I). Other explosion mod- 
els from that family show variations in the total 56 Ni 
mass, explosion energy, morphology and ejecta dynam- 
ics. With the addition of even more advanced physics, it 
is natural to expect that interesting new behaviors will 
emerge. Using a more realistic progenitor structure and 
ignition conditions may affect the nucleosynthesis in the 
model, thereby altering the composition of the core re- 
gion, energetics of the deflagration, and overall explosion 
energy. That in turn is likely to influence certain model 
observables, such as the absolute luminosity, the disper- 
sion in _B-band decline rates, and the model continuum 
polarization. 

A major advance of the multi-dimensional radiative 
transfer calculations presented in this paper is our abil- 
ity to quantify the dependence of model observables on 
the viewing angle. We identified several intriguing ori- 
entation effects in the Y12 model which suggest specific 
ways in which the asphericity of SNe la may contribute to 
their photometric and spectroscopic diversity. In partic- 
ular, asymmetry in the velocity distribution of iron group 
elements leads to a strong dependence of the i?-band de- 
cline rate on viewing angle. This fact both constrains the 
degree of compositional asymmetry present in the ejecta 
of SNe la and points to a potential source of dispersion 
in the SN la width-luminosity relation. Although not the 
case for the Y12 model, one may be tempted to wonder 
whether a multi-dimensional SN la model could be imag- 
ined such that the peak luminosity variations of SNe la 
and their correlation with the light curve width might 
arise, at least to some degree, from an orientation effect. 

The comparison of SN theory and observation faces 
new and interesting challenges once multi-dimensional 
models are considered. Given the viewing angle depen- 
dence, the predictions of aspherical models bear an in- 
trinsic multiplicity. Model validation can then no longer 
be limited to the traditional exercise of matching syn- 
thetic light curves and spectra to individual SN obser- 
vations. Rather, we must also study the probability 
distributions and dispersion levels characterizing various 
model observables (e.g., peak magnitudes, decline rates, 
line velocities and polarization levels) along with the in- 
ternal correlations relating different sets of such observ- 
ables. These statistical properties of the model can then 
be compared to those of a relevant sample of observed 
SNe. Although we have touched on this methodology in 
several places here, more generally the situation will be 
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further complicated by the presence of additional dimen- 
sions of diversity other than asymmetry alone. For ex- 
ample, the DFD and other theoretical paradigms admit 
a family of SN la models which range in 56 Ni mass. More 
extensive validation studies thus await further specificity 
of the theory. Equally important will be the development 
of massive and easily accessible observational databases 
containing light curves, spectra, and polarization for a 
large number of SNe. 

The general success of the Y12 model in reproducing 
the basic properties of observed SNe la is not entirely sur- 
prising, considering that the DFD scenario represents a 
variation on the standard ID delayed-detonation models 
which have bee n shown to be in fair agreement with ob - 
servations fe.g- lHoflich fc Khokhlo^[l^lHolfirchlll995f ). 
Moreover, the ejecta asymmetry studied here may be 
generic to a broader class of SN la explosion models. 
Standard delayed-detonation models, for example, show 
a similar systematic offset in the ejecta compositional 
structure if the trans ition to detonation is triggered a t 
one point off-center (|Livnd Il999l ; iGamezo et al.l [20051) . 
In such models, nuclear burning extends to higher veloc- 
ity on the side of the ejecta where the detonation takes 
place, just the opposite of that seen in the Y12 model. 
Nevertheless, the overall morphologies are similar, and 
the models likely share many of the same photomet- 
ric/spectroscopic variations and polarization signatures 
discussed here. 



Ultimately, the viability of the DFD explosion scenario 
is subject to further tests of the robustness of the DFD 
mechanism itself once the simplifying assumption of axial 
symmetry is removed and when more realistic initial igni- 
tion conditions are considered. Nevertheless, the overall 
good agreement of the Y12 model with a broad range 
of SN la observations certainly warrants future, more 
in-depth investigations into the DFD and related multi- 
dimensional delayed-detonation scenarios. This agree- 
ment is especially impressive given that the model was 
not artificially manipulated in order to match observa- 
tions, rather the system was evolved without intervention 
from the initial ignition to several months past maximum 
light. We believe that this sort of end-to-end methodol- 
ogy is critical in developing and validating future theo- 
retical models of thermonuclear supernova explosions. 
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